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' ROLL-UP OF A VORTEX SHEET* 

R. A. Mitcheltree** 
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Abstract 
■ ■ % 

The roll-up of a vortex sheet is analyzed by two approaches. The first 
is basti on the exact compressible Euler equations while the second is based 
on the exact incompressible Navier-Stokes equations. The inviscid calcula- 
tions for the two-dimensional problem do not indicate any roll-up of t! e 
sheet. On the other hand, the viscous calculations capture the dynamics 
of the roll-up rather well. This suggests that the generally held views 
regarding the roll-up process of aircraft wakes, namely, that it be treated 
as an inviscid process, may not be completely accurate. 

Introduction and Approach 

This research represents the results of the first phase of a research 
effoit dealing with aircraft wakes. An excellent review of this problem 
is given by Donaldson and Bilanin.^ According to Reference 1, understanding 
aircraft wakes entails the understanding of four problems. These are: roll- 

up, interaction and stability, aging and atmospheric effects. 
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The object of this investigation is the roll-up region. Its length 
scale is of the order of the aspect ratio over the lift coefficient times 
the span length. As indicated in Reference 1, processes in this region are 
essentially inviscid. Because of this, we felt that an approach based on 
the exact Euler equations should reproduce the experimentally observed 
features of this region. Such an approach has a distinct advantage over 
other approaches reviewed in Reference 1 because the generality of the 
equations permits first, detailed study of the relevant dynamics processes 
and second, observation of the downstream evolution of vorcical structures 
in their proper spatial and temporal relationships. 

Although w? are dealing with low Mach numbers, the compressible form 
of the Euler equations is employed. This is a crucial aspect of the form- 
ulation because the problem is dominated by boundary conditions. Figure 1. 

One does not know before hand whether inflow or outflow conditions should 
be imposed at a given boundary. Therefore, without employing a procedure 
that selects the proper boundary condition at a given instant and a given 
point one cannot hope to achieve an accurate solution. The proper boundary 

conditions for the compressible Euler equations, which are hyperbolic, are 
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determined from the method of characteristics and these have been imple- 
mented here. 

The calculation of the development of a voitex wake from a flat sheet. 
Figure 2, is considered. This problem, which is characterized as "formid- 
able" in Reference 1, is analyzed using the fourth-order Runge-Kutta method 
3 

of Jameson et al. Use of this method of solution brings up another crucial 
aspect of the formulation, namely, numerical damping. It is extremely 
important that the effects of numerical damping be reduced to a minimum in 
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order not to mask the relevant physics of the problem. There are three types 
of damping that are employed in Reference 3: second-order, fourth-order 

and enthalpy. Because the stagnation enthalpy is not constant for the 
problem under consideration, enthalpy damping cannot be used here. Second- 
order damping mimics viscous effects and, therefore, was not used. There- 
fore, only fourth-order damping was employed in this formulation. 

When one employs the compressible equations to study a phenomenon where 

the compressibility effects are small, the system of equations becomes stiff. 

The stiffness is a result of the disparity in the eigenvalues of the system. 

4 5 

Matrix preconditioning * may be employed to alleviate the problem. There- 
fore, a secondary objective of this investigation is to evaluate the utility 
of such schemes. 

Because of the inability of the Euler equations to predict the physics 
of the roll-up problem, the code that was developed in Reference 6 to study 
shear layers was employed. The code employs a vorticity-velocity formula- 
tion of the incompressible Navier-Stokes equations. Further details of the 
method are given in Reference 6. 

Results and Discussion 

The problem considered is that shown in Figure 2. At time t = 0, 
the velocity difference between the two streams occurs across one computa- 
tional cell with the result that, at t = 0, the vorticity is zero everywhere 
except at that cell. Because enthalpy damping could not be used to speed 
up the convergence to the steady state, the matrix preconditioning (M.P.) 
procedure of Turkel^ was employed. Figure 3 shows the convergence history 
of the solution based on the Euler equations in the absence and presence of 


matrix preconditioning,. It is ."seen from the figure that M.P. helps speed up 
the convergence. 

Figure 4 shows the vorticity contours that result from calculations 
based on the Euler equations with N indicating the number of iterations. 

The total number of iterations employed corresponds to the number of 
iterations needed to allow a disturbance moving with the lowest velocity 
to propagate across the computational domain. No roll-up is indicated 
in the figure. Figure 5 shows the vorticity contours based on the Navier- 
Stokes equations. It is seen that the results indicate roll-up. 

In conclusion, the above results suggest that the physics of vortex 
roll-up is best described by the Navier-Stokes equations. It appears that 
if we hope to understand the behavior of aircraft wakes then we must rely 
exclusively on the Navier-Stokes equations. 
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fin u re 4 |d vorticity evolution (inviscid) 
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figure 5a v vorticity evolution (viscous) 
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figure 5b vorticity evolution (viscous) 




